Singular current response from isolated impurities in d-wave superconductors 
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The current response of a d-wave superconductor containing a single impurity is calculated and 
shown to be singular in the low-temperature limit, leading in the case of strong scattering to a 
1/T term in the penetration depth A(T) similar to that induced by Andreev surface bound states. 
For a small number of such impurities, we argue this low-T upturn could be observable in cuprate 
superconductors. 
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Introduction. The quasiparticle excitations near the 
line nodes in the d-wave superconducting order parame- 
ter of the cuprate high-T c materials lead to well-known 
singularities in the current response of such systems, re- 
sulting in nontrivial magnetic field dependences of the 
London penetration depth in the Meissner state [1], and of 
the specific heat in the Abrikosov vortex state [2]. These 
excitations are also responsible for the marginal thermo- 
dynamic stability of the d-wave state itself, which has 
been related[3, 4] to the famous linear-T temperature 
dependence of the penetration depth A(T)[5]. 

Recently, it has been pointed out by Walter et al.[6] 
and Barash et al.[7], that, in addition to effects arising 
from extended quasiparticle states alone, there are 1/T 
contributions to the penetration depth below a crossover 
temperature T m o ~ \J§q/ AoT c , where £o is the coherence 
length and T c the critical temperature, due to zero-energy 
surface bound states. These states arise within a semi- 
classical picture when a quasiparticle reflected from the 
surface experiences a change of sign in the order param- 
eter; the upturns therefore only occur when the surface 
normal of the sample makes an angle close to ir/4 with 
the crystal axes and hence the antinodal directions of the 
d-wave order parameter. A A(T) ~ 1/T term in the pen- 
etration depth has in fact been reported by Carrington et 
al.[8], and attributed to the paramagnetic current carried 
by such states which dominates at low T. The field and 
interface angle dependence also appear to be in rough 
agreement with the predictions of Barash et al. [7] . Up- 
turns correlated with disorder introduced by Zn atoms 
observed by Bonn et al.[9] seem to vary from sample to 
sample, and are not currently understood. 

Most works in this area have treated disorder within 
effective-medium approximations which predict broaden- 
ing of low-energy quasiparticle states by a residual rate 
7, which is also roughly proportional to the residual den- 
sity of states at zero energy iV(0). Impurities are thus 
assumed to "smear out" the nodes of the gap, and there- 
fore inevitably cut off the singular behavior. For exam- 
ple, in the work of Barash et al.[7], no upturn in A(T) at 
low T is observed if 7 > T m Q. In the impurity-dominated 
regime, a quadratic temperature dependence is generally 
to be expected [11]. Impurity physics is therefore thought 
to compete[4] with nonlocal[10] and nonlinear[l] effects 



in the Meissner state, and A(T) changes from linear to 
quadratic below whichever scale is largest. 

There may be reasons, however, to doubt results 
obtained for the penetration depth using the self- 
consistent T- matrix approximation (SCTMA), which de- 
fines a translationally invariant effective medium for the 
impurity-averaged dirty d-wave superconductor. First, 
it has been shown to break down in two dimensions [12] 
and corrections due to weak localization and correlated 
order parameter response have been under intense in- 
vestigation recently. [13] Secondly, as with any effective 
medium theory, it must break down when impurities are 
sufficiently isolated. Finally, there are clear indications 
that single impurity bound states in unconventional su- 
perconductors are analogous in some respects to surface 
Andreev states. [14] Here we point out that under some 
circumstances isolated strong impurities can themselves 
make singular contributions to the penetration depth and 
other thermodynamic quantities. While we can present 
only estimates for the effect of a thermodynamically finite 
density of dilute impurities, we perform an exact calcula- 
tion for the current response of a d-wave superconductor 
in the presence of a single strong scatterer. 

Current response for single impurity. We consider a 
pure d wave superconductor described by Nambu propa- 
gator G k (cj) = {iuJ n T + CkT 3 + A k Ti)/D k , where T> k = 
+ £ k + A k , £ k = e k — jU is the 1-electron spec- 
trum measured relative to the Fermi level, and the Tj 
are Pauli matrices. The d-wave order parameter on the 
model cylindrical Fermi surface parametrized by angle 
4> is A k = Ao cos 2<f). In the system with one single 5- 
function impurity of strength Vq located at position R, 
the Green's function is exactly given by 

<2 kk ,M = G k M£ kk , + e l ( k - k ')- R G k (c)T(c)G k ,(c), (1) 

where T is the exact T-matrix for the impurity. We now 
ask what the change in the current response due to a 
single impurity is in this system. Note that we intend to 
study at first only the linear response to an external vec- 
tor potential A a la Kubo, but in the "unperturbed state" 
described by the exact one-particle eigenstates with the 
single impurity present. This is simply the usual nonlocal 
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expression for the current response in Coulomb gauge, 



J(r) 



J dr'K(r,r' 



)A(r'), 



(2) 



with^(r,r') = -c/(47rA 2 (T))5(r - r') + 5K(r, r'). Here 
we take A(T) to be the unperturbed London penetra- 
tion depth of the <i-wave superconductor, with pure be- 
havior A(T) ~ A ( l + (log2)T/A ) for T < A and 
A = yjmc 2 /Anne 2 . The Fourier transform 5K_(p, q) is 
now the change in response due to the impurity, which 
may be easily expressed in terms of the exact 1-impurity 
T-matrix using (1), 



5K^(p, q) = -"-e^-^ T J2 E \^ ( 3 ) 

u> 1 

{vf_ f vf_ 5 G?(o;)G?_ q ( £( ;)T(a;)G?_ p H 
+vf +? vf +§ G? +p (o;)T(a;)G? +q (a;)G?( W ) + 

k 

vr +t v^pG k (c,)T( w )G^ q ( W )G?( W )r(c,)G°_ p (c,)}, 

where w is an internal fermion Matsubara frequency, V is 
the volume of the system, and Vk = de/dk is the electron 
velocity. 

We now need to solve (2) together with Maxwell's 
equations to determine the spatial dependence of the vec- 
tor potential caused by the combined effects of Meissner 
screening and impurity scattering. Since we are primar- 
ily interested in long-wavelength effects, the perturba- 
tion may be considered to be of order 1/N, where TV 
is the number of atoms in the crystal; we may there- 
fore clearly treat the problem pcrturbativcly, by writing 
A(r) = A (z) + <5A(r), where the unperturbed solution is 
taken to be the London result (we assume Ao >> Co for 
the homogeneous system) A (z) = A (0) cxp — |z|/Ao- 
Note z > is the coordinate describing the theorist's 
half-space of superconducting material, and the solution 
A(r) is extended to unphysical values z < as an "image 
vector potential" to allow a solution by Fourier trans- 
form. Specular scattering of quasiparticles from the ac- 
tual surface is assumed throughout. 

To linear order in <5A, the problem can be cast as a 
differential rather than integro-differential equation, 



V 2 <5A 



1 



A(T) 



XA 47r 

-5A = — 



J dr'SK(r,r')A {z') (4) 



whose solution can be obtained by Fourier transform, 

F(q) 



F(q) = y dre l * r J dr'SK{r,r')A {z'). (5) 

With the expressions (2) and (5) it is clearly straight- 
forward, if tedious, to evaluate the full spatial depen- 
dence of the vector potential or the current around the 



impurity site. It is clear from physical considerations 
that the current at any given point in space r will de- 
pend both on the distance from the surface z and on the 
distance from the impurity r — R|. In general, then, the 
currents and field B = V x A are functions of x, y and z 
everywhere. We postpone the full evaluation to a subse- 
quent study. At present we would like simply to obtain 
an estimate of the size of this effect, which we do by ex- 
amining the perturbation of the component of the field 
along the initial applied field B(0) which we take to be 
along y. The penetration depth is usually defined, even 
in cases where the decay of the fields is not exponential, 
i.e. nonlocal electrodynamics, as Atot = J dzB y /B y (0). 
So for our current purposes we will simply adopt this 
definition and evaluate B y (x = y = 0,z), i.e. at the 
transverse position of the impurity R = (0, 0, z ). In this 
approximation we find 



5X SA x (0) + J^°dzd x 6A z + Xd z SA x {0)->dJA z (0) 



A 



A 0x {0) 



(6) 



A careful examination of these terms and glance at Eq. 
(5) shows that the size of <5(1/A 2 ) ~ 25X/X 3 is simply 
set by SK_(p, q) with p,q ~ 1/A, as intuitively expected. 
Expanding Eq. (3) for p,q <C kp and performing the 
integration over energy £ k yields from the first two terms 



T f 2 * 2 W A 2 (12£ 2 + 7? 2 ) 



V 



(7) 



where & = ^Juj 2 + A 2 , T - (^ )- 1 G /(c 2 - G 2 ), with 
Go - (l/27rJV )Tr£ k G_ k = (2iu/^ 1 N )K{A /^), is 
the To component of the T-matrix, and c is the cotangent 
of the s-wave scattering phase shift (c = corresponds 
to infinitely strong scattering). K, is the complete elliptic 
integral of the first kind. Note that the third term of 
(3) contains two factors of Tq and therefore formally of 
order 1/N 2 ; we therefore neglect it in its effect on long- 
wavelength properties like <5A, although it will contribute 
to local properties. The quantity i] p = v-^ ■ p ~ vp/X, 
where the last approximation follows because the primary 
spatial gradients are perpendicular to the surface with 
normal z, and because we are interested in typical p ~ 
1/A. In Figure 1, we show a numerical evaluation of 
(7) for three different values of the impurity scattering 
phase shift. As the scattering approaches the unitarity 
limit, the perturbation to the penetration depth is seen 
to diverge with decreasing temperature. 

To get an analytical estimate of the upturn in the uni- 
tarity limit, we first perform a nodal expansion of the 
order parameter, A k — 2Ao(</> — </>„), where </>„ = tt/4, 
37r/4, leading for T <C A to the result 



JVn" 



1 



4ttA 2 27r 2 TVn log^- 



(8) 



where the infrared divergence in (7) was cut off by n p 
or cA , and no is the appropriate minimum Matsub- 
ara index no = max(?7p/7rT, cAq/ttT, 1). For a (001) 
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FIG. 1: Normalized change in penetration depth 
8\/[\o{Nq /VT c )] due to single isolated impurity vs. 
normalized temperature T/T c for three values of the cotan- 
gent of the impurity phase shift c = 0,0.1,0.3. For this plot 
we took ri P ~ 0.01T C . 



surface in the cuprates, r/ p ~ (v F /vp~)(£,o/Xo)Ao — 
10~ 3 — 10~ 4 , so we get a divergence which goes as 
5X ~ X a E F /[NT\og(A Q /T)}. The factor of 1/N arises 
since (8) is proportional to N^ 1 /V ~ a 3 /V, where a is 
the lattice spacing. The estimate holds only for temper- 
atures down to T* = max(?7p,cAo) at which the diver- 
gence is cut off, of order tenths of a Kelvin or less in the 
cuprates if c is taken to be 0. 

Finite density of scatterers. We would now like to ex- 
tend the above estimate for a single impurity to the case 
when a finite density of impurities is present. As men- 
tioned above, when brought into proximity these impuri- 
ties interfere with one another via hybridization of quasi- 
particle bound state wave functions, leading eventually 
to an "impurity band" and residual density of states at 
the Fermi level, at least in 3D [13]. In this situation the 
SCTMA is the appropriate approximation, and a T 2 be- 
havior in the penetration depth is found rather than an 
upturn. [11] The conditions for the formation of this band 
are subtle, and we do not address these questions here. 

Intuitively, it seems likely that impurities will act as 
independent sources of current distortion if they are suf- 
ficiently far from one another, but it is not so sim- 
ple to specify what "sufficiently far" means. We note 
that both the distortion of the current distribution aris- 
ing from 5K(r,r') and the impurity bound state wave 
functions gcncrically decay over lengths of a few £ - If 

— 1/2 

the typical interimpurity spacing I = n i 1 a (in 2D) 
is much greater than £o; each impurity in a penetra- 
tion length A contributes 8X to the observed <5A tot . On 
the other hand, even if £ « £ , there is a probability 
P(R) ~ exp— R 2 /2£ 2 that an impurity will find itself 
with no other impurities within a radius R, if all impuri- 
ties are randomly distributed. We therefore expect that 
an upper bound to the observable penetration depth up- 
turn in the case of randomly distributed unitarity limit 
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FIG. 2: Change in penetration depth <5Atot(A) due to isolated 
unitarity limit scatterers vs. temperature T(K) for three im- 
purity concentrations n t = 10~ 2 , 10~ 3 , 10 -4 . A /T c = 2.14, 
Vp/Tc = 0.001, T c = WOK. Dashed line: T-variation of 
pure penetration depth (log2)T/A . Inset: 8X tot vs. n t at 
T = IK. 



scatterers will be 



SX 



tot 



47r 2 flog(2A /7rT) 



(9) 



with T = T + T* . The estimate is an upper bound be- 
cause we do not currently know how to calculate accu- 
rately a geometric factor arising from the contribution 
of those impurities which arc located such that the tails 
of nodal quasiparticle resonant states (which do not de- 
cay as e - ^°/ r but rather as 1/r at lu = 0)[15, 16] over- 
lap. This factor will lead to significant hybridization and 
smearing of the singularity at impurity densities below 
that given by the criterion above. We nevertheless sub- 
stitute some numbers into (9) to check for plausibility, 
taking E F /A ~ 10, A = 1500A, and £ = 25A for the 
cuprates. In Figure 2, we plot the results of (9) in labora- 
tory units given these assumptions. It is seen that, within 
this crude approximation, the size of the upturn expected 
is not monotonic with impurity concentration, but at a 
temperature of IK peaks around rij ~ 10~ 3 . The position 
T 2 * of the minimum in A tot (T) will now be set by com- 
paring the impurity-induced change in A with that due to 
thermally excited quasiparticles, <5A(ni,T 2 *) ~ T 2 *-^o/Ao; 
for impurity concentrations of rij = 10~ 3 — 10~ 2 , this 
temperature will be of order 1-2K, as seen in Figure 2. 

The size of the increases predicted seem quite reason- 
able, and are of the same order of magnitude as the up- 
turns seen in the work of Bonn et al.[9] and Carrington 
et al.[8] Although the latter authors observed a subtan- 
tial increase in the size of the signal when the proportion 
of (110) surface in the given sample was maximized, in- 
dicative of an Andreev surface bound state contribution, 
there was a substantial signal even in samples with only 
(100) surfaces, where surface bound states should not ex- 
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ist. Although the authors' suggestion that this result is 
due to (110) faceting in these samples is possible, the 
mechanism we suggest may also be present. 

While the estimate (9) is based on the assumption of 
randomly distributed impurities, it is important to recog- 
nize that in real systems clustering will take place. In this 
case even a relatively disordered system may have signif- 
icant numbers of isolated impurities or isolated atomic- 
scale clusters which give rise to a singular current re- 
sponse. In this case we might expect significant sample- 
to-sample variation in SX in samples with identical aver- 
age concentration rij. This appears to be precisely the 
effect found by Bonn et al. in their measurements on 
YBCO single crystals with 0.31% Zn[9]. 

Magnetic field. Since the penetration depth upturn in 
the case considered here is due to the large number of 
quasiparticle excitations near the nodal directions, as is 
the similar upturn in the case of Andreev surface states, 
we might expect any physical effect which smears the gap 
nodes to cut off the upturn. In particular, the orbital 
coupling to an applied magnetic field (nonlinear electro- 
dynamics) will suppress the upturn as it does in the An- 
dreev case [7]. We need only add the Doppler shift of 
the quasiparticle energy in any of the expressions above, 
iuj — > iu> + v s • k, where v s is the local superfluid ve- 
locity. A typical shift v s Uf — (H/H )A , where H is 
the applied field and H = 3<&o/(7r 2 £ Ao) is of the or- 
der of the thermodynamic critical field, and the effect of 
the field may be included approximately by generalizing 
f -> T + Tj* + (H/H )A in (9). We postpone a more 
quantitative study of the field dependence to a later work. 

Conclusions. We predict that a ci-wave superconduc- 



tor with isolated strong nonmagnetic impurities should 
exhibit an upturn in the penetration depth varying as 
1/TlogT below a temperature scale set by the disorder 
in the system and above one set by the bulk penetra- 
tion depth or the scattering phase shift. This effect may 
have contributed to upturns already observed in exper- 
iments where Andreev surface states, which produce a 
similar effect, are not indicated[8, 9]. The physical origin 
of the upturn in both cases is the depletion of the su- 
percurrent at low T in low-energy bound states, but the 
impurity-induced upturn is independent of surface geom- 
etry. We argue that the magnitude of the effect should 
have a characteristic dependence not only on the con- 
centration of impurities, but on the exact nature of the 
statistical distribution of impurities in the sample. We 
emphasize that the estimates presented here are merely 
a crude plausibility argument that if such upturns are 
observed at low temperatures, they could arise from sin- 
gular current response of local defects in the crystal, and 
that it will be useful to study the dependence on disorder. 
Our work also allows in principle an exact calculation of 
the local current flow around an impurity; it will be in- 
teresting to see how the singular quasiparticle currents 
are distributed in the neighborhood of the impurity site, 
and whether these features can be detected by STM. 
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